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We report first principle numerical study of domain wall (DW) depinning in two-dimensional 
magnetic film, which is modeled by 2D random-field Ising system with the dipole-dipole interaction. 

We observe nonconventional activation-type motion of DW and reveal its fractal structure of DW 
near the depinning transition. We determine scaling functions describing critical dynamics near 
the transition and obtain universal exponents establishing connection between thermal softening of 
pinning potential and critical dynamics. We observe that tuning the strength of the dipole-dipole 
interaction switches DW dynamics between two different universality classes corresponding to two 
distinct dynamic regimes, motion in the random potential and that in the random force. 


Introduction - Motion of domain walls in magnetic 
nanowires and films is a key component of operation 
of any magnetic memory and logic device |l|-|3(. To 
a great extent DW dynamics is governed by pinning- 
depinning processes which control the operational speed 
and power consumption of a device and thus play cen¬ 
tral role in device performance [3, H|. There has been 
remarkable progress in description of pinned DW dy¬ 
namics based mostly on the elastic manifold model in 
a random environment [§]. A key property of such a sys¬ 
tem is that at zero temperature it experiences the dy¬ 
namic phase transition (depinning transition): At small 
external drives, F ^ F c , where F c is the critical pin¬ 
ning force, DW is immobilized (pinned) by disorder, 
whereas at F > F c it acquires a finite velocity v. The 
threshold depinning force F c is a critical point in a sense 
that at F > F c , the velocity exhibits critical behavior 
v ~ ( F—F C )P , as was proposed by Fisher[3 in the context 
of depinning of charge density waves. At finite temper¬ 
atures the velocity is always finite, and at F <C F c the 
domain wall exhibits highly nonlinear glassy response, 
so-called creep dynamics, with v oc exp(— const/F M )[(f]. 
The depinning transition gets rounded and acquires a 
meaning of the intermediate region separating the low 
force creep dynamics from the asymptotic linear response 
v oc F at F F c . The critical depinning behavior has 
to include temperature dependence and was conjectured 
to be of the form: v ~ 'F [(F/F c — 1) /T v ] @, 11]. 

All the above results were obtained within the elas¬ 
tic manifold approach, where the DW was modeled as 
an elastic membrane (in 3D) or an elastic string, if we 
discuss one-dimensional DW in a magnetic film. How¬ 
ever fundamental and successful, this description misses 
important processes that may become essential for the 
DW dynamics at elevated temperatures. Domain wall is 
multi-valued, so ahead boundaries can merge with the 
‘main’ interface behind. 

This poses a challenge of developing first principle ap¬ 
proach starting from the microscopic model that captures 
basic physics of the magnetic system. Taking up upon 


this challenge we reveal the fractal structure of the do¬ 
main wall and uncover the critical dynamics at the depin¬ 
ning transition, and relate the observed critical exponents 
to those of finite temperature creep dynamics. Further¬ 
more, we uncover the role of strength of the dipole-dipole 
interaction in determining the proper dynamic universal¬ 
ity class. 

Model and method - We model the two-dimensional (2D) 
magnet subject to quenched disorder by the 2D random- 
field Ising model with the dipole-dipole interaction, and 
the dynamics is controlled by the external driving field: 

H = -jY J SiS j + V dd Y J ^4 l -Y.( hi+H ) Si ’ (!) 

(i,j) i<j r U i 

with Si = ±1 at site i. The first term of the Hamilto¬ 
nian is the ferromagnetic coupling between one spin and 
its nearest neighbors. Hereafter we measure the energy 
in the units of the coupling J. The second term is the 
magnetic dipole-dipole interaction with nj = \i — j\ and 
Vdd a parameter for interaction strength. The on-site 
random field hi distributes uniformly within an interval 
[—A, A] which generates random pinning potentials. H 
is a uniform magnetic field which drives the domain wall. 

Our simulations are performed on L x L square lat¬ 
tice. A flat domain wall between spin +1 and spin —1 
is created along y axis at x = 1 as the initial condition. 
The magnetic field is applied to drive the domain wall 
in the positive direction of x axis in accordance to the 
Metropolis algorithm with single-spin flip EH- Periodic 
boundary condition (PBC) is adopted at the domain-wall 
(y) direction, whereas Anti-periodic boundary condition 
(APBC)[I3 the moving (x) direction. The number 
of independent runs is at least 3000. The time unit is 
defined by a sweep of Monte Carlo trials over the whole 
system, and the velocity is defined by v = dM/2Ldt in 
steady states, with M the total magnetization. For a 
small system under large driving field, it reaches a steady 
state quickly [t ~ 10 for warm-up and ~ 10 2 for statis- 
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so do the corresponding values of the depinning field. 

It is observed that for a given field and system size, the 
domain wall may either be pinned inside the sample, or 
it may go through from edge to edge. We call the latter 
case a depinning event and evaluate the corresponding 
depinning probability Pq. To determine a true value of 
H Cl one thus has to perform the finite-size scaling analysis 
of Pq which would contain H c as a parameter. As shown 
in Fig.QJ,, Po(L,H) increases sharply as function of the 
magnetic field in the interval H = 1.1 ~ 1.3. The curves 
corresponding to different system sizes cross at point of 
Po = 0.38 ± 0.04 at H = 1.214 ± 0.006. This determines 
the depinning field which does not depend on the system 
size and thus can be taken as a depinning field H c of a 
macroscopic system. 

The problem of depinning at zero temperature is in¬ 
timately related to percolation problem [b|. We thus 
assume that the depinning probability function has the 
form characteristic to the percolation problem (lij: 

P 0 (L,H) = *[(H/H c -l)L 1 /% (2) 


FIG. 1. (Color on line), (a): The probability Po versus the 
driving field with different system sizes and the corresponding 
scaling plot (inset). Vdd/J = 0.1 and A = 1.5J are used 
all through this work, (b): v — H characteristics at zero 
temperature with difference system sizes. The solid line is 
the fitting function (Eq. l|3|)) for L = 1024. 


tics), whereas typically 10 3 time steps for warm-up and 
statistics with regard to large systems under critical driv¬ 
ing field. The time scale in MC technique should be pro¬ 
portional to the real time, but a straightforward relation 
is not easy. In order to derive the correspondence, one 
need to compare the simulation results and experiments 
at least once. 

The system adopted in the present work is a coarse¬ 
grained one. The on-site spin is represented in terms of 
a block spin which contains n z x n xy x n xy unit cells 
in a thin magnetic film, with n z and n xy the number 
of unit cells along the out-of-plane and in-plane direc¬ 
tion, respectively. Then the energy unit Ej = n z n xy aA , 
with a the lattice constant of a real material and A 
the corresponding exchange stiffness. For Nd 2 Fei 4 B, 
a = 0.88 nm and A = 7.7 pJ/m[ll|, whereas a = 0.25 
nm and A = 10.3 pJ/m for Co layer in Pt/Co/Pt thin 
films fill. fl2j . The dimensionless temperature T in the 
present work can be related to the real temperature 
through: T ■ Ej/ks , with fee the Boltzmann constant. 
Taking n z = n xy = 1, one can have T = 0.1 approximates 
to 49 K for Nd 2 Fei 4 B, whereas 19 K for Pt/Co/Pt thin 
film. 

Zero-temperature depinning - To come up with the quan¬ 
titative description of depinning, we have to know its key 
characteristic, the zero-temperature depinning field H c . 
Finding its true value is a challenge since in finite sys¬ 
tems realizations of the random potential fluctuate, and 


where v is a universal exponent. By choosing the vari¬ 
able (H/H c — 1 )L 1 / V and v as an adjustment parameter, 
we find that at v = 1.33 ± 0.05, all the data points of 
P 0 (L,H) collapse onto a single curve as shown in the 
inset of Fig.[ljr. This procedure defines the exponent v. 

Now we are equipped to study the v-H characteris¬ 
tics in the depinning regime. We start with the zero- 
temperature behavior. The results for v-H for systems 
of different sizes are displayed in Fig.QJ). One sees that 
for L ^ 128 the curves do not practically depend on the 
size of the system. Assuming the standard v-H depinning 
relation [l5f 

v = v 0 {H/H c ~lf, (3) 

where /3 is a universal exponent, one finds vq = 1.16 ± 
0.01, H c = 1.214T0.006 and /? = 0.36T0.01 for L = 1024. 
This value is in a fair agreement with the 0 = 0.31 result 
obtained in two-loops RG calculations [RJ showing that 
the elastic manifold approximation works pretty well at 
zero temperature. 

Finite-temperature depinning - Now we turn to our main 
task, the finite temperature motion. To reduce the com¬ 
putation time we choose L = 512 system. Figured 
shows the expected increase in velocity at the given field 
upon increasing temperature and an appreciable tail be¬ 
low the depinning field due to thermally activation pro¬ 
cesses. 

We use the standard scaling ansatz fl7i - [l9j : 


;(T, H) = T 1 / <5 'F \(H/H c - 1 )T~ 


1//3S 


(4) 


with T(x) ~ x 13 as x —> oo. We achieve the best collapse 
of the data to a single curve with <5 = 2.76 ± 0.02 by 
adopting the values of H c and /? determined above, see 
Fig. [2]). Note, that these results cease to hold for large dd 
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FIG. 2. (Color on line), (a): Finite temperature v — H characteristics for Vdd/J = 0.1. (b): Scaling plot of v — H curves 
as vT~ x ^ vs. (H/H c — 1 The inset shows the same scaling behavior replotted in the semi-log scale, (c): Finite 

temperature v — H characteristics for Vdd = 0. (d): The corresponding scaling plot for Vdd = 0, the inset shows the same data 
in semi-log scale. 


interaction, Vdd/J > 0.5 where the ferromagnetic order 
is broken. 

At x < 0 the scaling function exhibits the asymptotic 
behavior d>(a;) = 1.05 exp [—0.70(—a;) 5 / 3 ], see Fig.[2j), 
and one arrives at the dynamics of the domain wall across 
the transition [13] given by: 


v = VxT 1 /* 




(5) 


where E c ss 0.81 is an energy barrier which governs the 
domain-wall velocity at finite temperatures, and the con¬ 
dition (35 = 1 is taken into account. Notably, the domain- 
wall motion is not the conventional Arrhenius-type. The 
origin of this nontrivial temperature dependence is the 
renormalization of the random potential landscape by 
thermal fluctuations. 

The exponent 5/3 is in an excellent ag reement with 
the prediction of Ref. @ (see Ref. [2(| for thermal 
depinning of vortex systems). This establishes an inti¬ 
mate connection between the critical depinning behav¬ 
ior of the domain wall at finite temperatures and ther¬ 
mal softening of the pinning potential. Furthermore, 
juxtaposing our v-H curves with those obtained ear¬ 
lier for the case Vdd = 0[l7j], one makes a remark¬ 
able observation. Figure^ shows the v{H) dependen¬ 


cies that at the first glance are not that different from 
those of Fig.[2ja. The scaling treatment, however, yields 
v oc exp [—U c { 1 — H/H c ) 5 / 3 /T ], see Fig.[2Ji, i.e. the Ar¬ 
rhenius activation behaviour with the barrier that scales 
as U c ( 1 — H/H c ) 5 / 3 , where U c ~ 0.67 is the bare energy 
barrier, H c = 1.289, (3 = 0.33, and (36 = 5/3. We now 
recall that this kind of thermally activated behavior at 
Vdd = 0 is expected for the case of the motion of DW 
in the field of the random force0, where an impurity 
“knows” to which magnetic domain it belongs in. This 
implies that the pinning energy barrier comprises the en¬ 
ergies of all pinning sites located in the area spanned by 
the domain wall during an elemental activation jump. 
In this case the contribution of the fluctuation thermal 
broadening of the domain wall position gives the negligi¬ 
ble contribution into the total pinning energy thus ther¬ 
mal fluctuations cannot appreciably reduce the depin¬ 
ning field. Hence the thermally activated motion retains 
its Arrhenius-like character. We conclude that varying 
Vdd one tunes the system between the random force and 
random potential pinning behaviors. Hence our findings 
provide an irreplaceable tool for identifying these distinct 
pinning mechanisms in the experiment. The behaviors 
summarized in Fig. [2] constitute the main results of our 
work. 

Another comment in order is that at finite temper- 
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FIG. 3. (Color on line), (a): Time evolution of domain wall near the depinning threshold at zero temperature, (b): Height- 
difference correlation function C(r,t) versus r for different values of t near the depinning threshold, (c): Saturated value C s (t) 
versus t. (d): Scaling plot for data in (b) with the dashed line y ~ x 2 ~ df . (e): Log-log plot of N versus l/l at different times, 
with l the ruler size and N the measured length in units of l. Parameters used in calculations are: L = 1024, Vdd/ J = 0.1 and 
H = 1.22. 


atures the definition of the depinning field H c is not 
straightforward. Our approach offers a systematic way 
for analyzing data at finite temperatures yielding H c and 
the thermal activation energy barrier E c simultaneously. 
Importantly, E c depends not only on the strength of ran¬ 
domness A but also on the competition of the exchange 
coupling J and the dipole-dipole interaction strength Vdd- 

Domain-wall morphology - Next we investigate the DW 
morphology during the depinning process for the system 
with Vdd/J = 0.1. To this end we set a flat domain 
wall along y axis at x = 1 at t = 0 with H = 1.22 
just above Hq = 1.214, and drive it along x direction at 
zero temperature. As shown in Fig.[3Ji, the domain wall 
evolves rougher with time, and develops a fractal struc¬ 
ture. Moreover, there remain several small unflipped- 
spin areas (black puddles) forming the “ lakes” inside 
the domain. The multiconnected nature of the flipped 
domain originates from spatial fluctuations of the depin¬ 
ning field due to random character of the pinning poten¬ 
tial: there are lacoons where the local depinning field still 
exceeds the driving field. The “overhangs” of the frontier 
of the domain wall are of the same origin. To quantify 
this multi-valued domain-wall morphology, we define the 


modified height function h(y, t ) [22j 

L 

%>*) = ( 6 ) 

X — l 

with 9{x) the unit step function and S x , y (t) the spin value 
on the site (x,y) at time t. The function h(y,t) describes 
the total number of flipped spins along the line y at time 
t. It is obvious that h(y,t ) describes the position of do¬ 
main wall if there is no “ lakes” and “ overhangs”. 

The height-difference correlation function C (r, t) de¬ 
scribes the domain-wall roughness characteristics and is 
define asjH, [Hj]: 

C(r,t) = ^([h(y + r,t)-h(y,t)] 2 y (7) 

with r the distance between two points in y direction. 

As shown in Fig. [3}y for fixed t , C(r) increases from 
zero with r and saturates at large r, and the saturated 
value C s (t) increases with t. These properties can be un¬ 
derstood from the time-evolution process of domain-wall 
morphology as displayed in Fig.[3J.i. The initial domain 
wall is a straight blue line with no height difference. By 
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applying driving field, locally meandering segments ap¬ 
pear first (see the light blue regions in Fig.[3Jj). Correla¬ 
tion of domain-wall positions only exists in a small length 
scale. As time evolves, the meandering segments spread 
out along both the domain-wall (||) and moving (_L) di¬ 
rections, leading to rougher structures. There are two 
f-dependent correlation lengths: £\\{t) and £j_(t), grow¬ 
ing with time as £_l(£) ~ £||M a ~ with z the dy¬ 
namic exponent and a the roughness exponent, and the 
correlation function evolves as|22l|: 

0(r,t) ~ UMr/Z || W] - (8) 

with g{x) saturates at constant as a; >> 1. 

As displayed in Fig. [3fc, we obtain a/z = 0.82 ± 0.01 
in terms of C s (t) ~ f“/ z in the large r limit of Eq. <(8|). 
By choosing z = 1.43 ± 0.01, all the data collapse into 
a single curve as displayed in Fig.[3ji, which determines 
the dynamics exponent 2 . The roughness exponent a is 
then estimated as a = 1.17 ± 0.02. 

We then study the fractal geometry of the domain 
wall. For a fractal structure, the measured length N 
in units of ruler size is related to the ruler size l by: 
df = log N/ log(l/7), with df the fractal index. Through 
log-log plot of N versus l/l as shown in Fig.[3}3, we ob¬ 
tain df = 1.25 ± 0.01. We notice that C(r,t) ~ r D ~ d f 
only holds for r < £||(£)[22|. As shown by the dashed 
line in Fig. Oil, the exponent 2 — df ~ 0.75 appears in the 
small scaling variable r/t 1 11 ^ limit of function g{r/t l / z ). 
Importantly, the fractal morphology is a signature of the 
depinning region. At large drives H, the DW gets flat [23] . 

Discussions and conclusions - To conclude, we have in¬ 


vestigated the depinning dynamics of magnetic domain 
wall with the dipole-dipole interaction, and obtained five 
critical exponents. It is confirmed that the scaling rela¬ 
tion fi/v = z — a[24|] is satisfied. The case without the 
dipole-dipole interaction has been investigated previo usly 
by MC simulationsfl?) and Langevin dynamics [25- 26|. 
The critical exponents are different for the two cases, 
indicating the two cases belong to different universality 
classes. 

The existence of two universal classes in elastic man¬ 
ifolds in random potentials has been addressed by two 
of the present authors |To|| for vortex dynamics in type- 
II superconductors, where Bragg glass and amorphous 
vortex glass (AVG) correspond to weak and strong ran¬ 
dom pinning potentials. In heavily disordered AVG the 
depinning dynamics is of Arrhenius-type, while in the or¬ 
dered Bragg glass state one observes non-Arrhenius-type 
behavior, similar to the present domain-wall system. In¬ 
terestingly, scaling functions of two universality classes in 
vortex dynamics share a simple exponential form, which 
differs from our results for the different dimensionality of 
the space and the elastic manifold; this intriguing issue 
calls for further investigation. 
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